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Abstract 

We show that parametric level correlations in random-matrix theories 
are closely related to a breaking of the symmetry between the advanced 
and the retarded Green's functions. The form of the parametric level cor- 
relation function is the same as for the disordered case considered earlier 
by Simons and Altshuler and is given by the graded trace of the commuta- 
tor of the saddle-point solution with the particular matrix that describes 
the symmetry breaking in the actual case of interest. The strength factor 
differs from the case of disorder. It is determined solely by the Goldstone 
mode. It is essentially given by the number of levels that are strongly 
mixed as the external parameter changes. The factor can easily be esti- 
mated in applications. 



1 Introduction 

Parametric level correlations in chaotic and disordered systems have re- 
ceived much attention in the early 1990's (see the review ^ and references 
therein) . This development culminated in the seminal work of Simons and 
Altshuler |5| who showed that such correlations have a universal form and 
who calculated some of the correlation functions for disordered systems 
explicitly. 

In the present paper, I take a fresh look at this problem. This is mo- 
tivated by two circumstances, (i) The work of Simons and Altshuler does 
not address level correlations of random matrices (but rather of chaotic 
and/or disordered systems). But level correlations of random matrices 
do play a role in some applications of random-matrix theory. A case in 
point concerns correlations of levels with different spins in atomic and 
nuclear systems While the form of the correlation function obtained 
by Simons and Altshuler is expected to be unchanged, it is necessary to 
determine the dimensionless parameter which governs its behavior, and to 
connect that parameter with physical parameters of the system at hand. 
In particular, the concept of "level velocities" introduced by Simons and 
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Altshuler needs to be reconsidered. It will be shown that in contrast to the 
case of disorder, the strength parameter in random-matrix theory is not 
influenced by a coupling of the Goldstone mode with the massive modes, 
(ii) Parametric level correlations can be seen as a manifestation of sym- 
metry breaking. The broken symmetry is that between the advanced and 
the retarded Green's functions. I aim at a presentation which displays 
this fact as clearly as possible. With this insight, writing down the form 
of the correlation funtion is quite straightforward. 

The correlation functions will be given for the GOE and for the GUE. 
The results are also compared with the two-point correlation function for 
the GOE GUE transition caused by time-reversal symmetry breaking. 
We shall see that in the latter case, symmetry breaking acts differently. 



2 Formulation of the Problem 

The ensemble of Hamiltonians H has the form 

H = Hxcos{X) + H2S,in{X) (1) 

where X is & dimensionless parameter and where H\ and H2 are uncor- 
related random matrices belonging to the same symmetry class of one of 
Dyson's three canonical ensembles. We wish to calculate the parametric 
correlation function 



fc = trf— r ltr[ 1. (2) 

^E+-H{Xy ^E--H[X'y 

The overbar denotes the ensemble average. The function k contains quan- 
titative information about the way in which the spectra at parameter val- 
ues X and X' are correlated. For X = X' , k coincides with the standard 
two-point correlation function. 

Both in the case of disordered systems and in the present case, one 
needs to calculate k only for small values of |X — i.e., perturbatively. 
The reason is that we are interested in local (rather than global) changes 
of the spectrum. The former involve an energy scale of order d, the mean 
level spacing, the latter, an energy scale of order Nd where A'' ^ cxj is the 
dimension of the matrices Hi and H2- Then, the function k = k{e, X —X') 
depends only upon the difference e — Ei — E2 of the energies of the two 
Green's functions. 

I expand the Hamiltonians H(X) and H(X') in Eq. Q around the 
mid-point Xo = {1/2){X + X') in powers oi X - Xo = (1/2)(X - X') 
and of X' ~ Xq — (1/2)(X' — X), respectively, and keep only terms up to 
first order in {X - X'). Then, 

H{X) ^ Ho + {l/2){X -X')V , 
H{X') ^ Ho - {1/2){X - X')V , (3) 

where Ho = H{Xo) and where 

V = H2 cos(Xo) - Hi sin(Xo) . (4) 
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The random matrices H{Xo) and V are uncorrelated, 



HoV = . (5) 

This follows from the fact that Hi and H2 are uncorrelated, H1H2 = 0. 

To identify the small parameter of the expansion, I define (as usual) 
the spreading width due to the perturbation as 

= 2tv{X - X' fV^/d . (6) 

The spreading width is a measure of the energy interval within which the 
levels of Ho get strongly mixed as the parameter changes from X to X' . 
We are interested in values of which are of the order of d (rather than 
Nd). We normalize the variances of Hi and H2 in the usual manner. 



m^,AH,)^^ = - ; j = l,2-f,^u , (7) 

where fi and v are level indices, and where 2A is the radius of the semicircle. 
The mean level spacing of Ho in the centre of the semicircle is given by 
d — ttX/N. With these conventions, we have 

r^^2{x-x'fx. (8) 

We shall see that the dimensionles parameter which governs the level cor- 
relation function is given by r^/d = (2/n)N{X-X'f. For this parameter 
to be of order unity, we must have that {X — X')^ is of order 1/A^. This 
jusitifies our perturbation expansion and the fact that we keep only linear 
terms in {X - X'). 

Substituting H(X) and H{X') from Eqs. igj into Eq. © yields 

k{e,X ~X') 



- Ho - (1/2)(X - X')V^ ^^E^ -'H0 + (1/2)(X - X>)V^ ' 

(9) 

Eq. © displays explicitly the fact that the perturbation V breaks the 
symmetry between the retarded and the advanced Green's functions. This 
is essential for the supersymmetry calculation of k{e,X — X'). 



3 Supersymmetry 

The supersymmetry method has become a standard tool in random- 
matrix theory. Therefore, I confine myself to giving the essential steps in 
the calculation. I do so for the case where both Hi and H2 belong to the 
GOE, and give only results for the GUE. 

I proceed as in Ref. also use their notation, and arrive at the 
following form of the generating function, 

ZiEi,E2;X,X',J)= /d[*]exp{£(*, J)} , (10) 
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where the Lagrangian is given by 

£ = (l/2)i(*^L'/^L»''L'''^*) . (11) 

Here is a graded matrix of dimension 8, given by 

D'^ = (E- H + i(5 + J - (l/2)£:) . (12) 

According to Eq. J^J, the matrix H has the form 

H = J/(Xo)l8 + (l/2)(X-X')l/L . (13) 

Here Ig denotes the unit matrix in eight dimensions, while 

L = diag(l, 1,1, 1,-1, -1,-1,-1) (14) 

is the matrix which breaks the symmetry between the advanced and the 
retarded Green's functions. We want to calculate the two-point function 
and accordingly put 

J = 5M^diag(-jl,-jl,+jl,+jl,-j2,-j2,+j2,+j2) = 5M^(jl,j2) • (15) 

The last equation defines (ji,j2)- 

The ensemble average is given in terms of the second moment of the 
term {i/2){-^'^ V-'^iiL^^'^<^), 



[(i/2)(vE'tLi/2HLi/2^,)]2 

= -(AV(2iV)) Y f*L(i'/').c(L'/')c.*.. 



-(AV(8iV))(X-X')' J2 (*Ic«(i'/')««i««(i'^')cc*.. 



X [^ip{L''^)PPLp^{L'/^)pp<i>^^ . (16) 

The summation over a, 13 runs from 1 to 8, that over /i, v from 1 to A^. I 
define 

+ (l/8)(A-A')'iA^L,,(L'''')„cVM-V'i/3(i''")/3/3i/3/3 • 

/J 

(17) 

This equation clearly displays the separate contributions from Ho and 
from the symmetry-breaking term VL. Under neglect of higher-order 
terms in {X — X')^ (which we have shown to be negligible for N oo), 
the right-hand side of Eq. 1)16^ can be expressed in terms of A, yielding 



[(i/2)(*tLi/2HLi/2vI/)]2 = ^trgj^^) . (18) 
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The Hubbard-Stratonovitch transformation yields now for Z the form 

Z(i5i,i52;X,X',J) = y d[a]exp|~^trgJf72)-^trg,lnN(J)| , 

(19) 

where 

N(J) = £18 - (l/2)£: + z5 - AE + (ji, ja) (20) 

and 

E = (7+ (1/8) (X - X')^ icrL . (21) 

We use the saddle-point approximation, omitting terms which are of order 
1/A^. These are the terms proportional to £, to [X' — X)^, and the source 
terms. The saddle-point equation 



A 



E\ii - Act 

has the standard solution 



(22) 



aa = Tg V^To (23) 



with (jj) diagonal and given by 



a?, = ^ - iAoL (24) 



and Ao = a/1 - (_B/(2A))2. The full sigma matrix is written as 

a^aG + Sa = aG + T^^SPTo . (25) 

It remains to work out the integrals over the massive modes, and over the 
Goldstone mode. 



4 Integration over the Massive Modes 

In applications of the supersymmetry formalism, one would normally skip 
the present Section because the integration over the massive modes is 
known to simply yield a constant. However, in the work of Simons and 
Altshuler |2, it is shown that the strength of the parametric level cor- 
relation function depends upon contributions due to the coupling of the 
Goldstone mode with the massive modes. Is such a mechanism also opera- 
tive in the present case? To answer this question, I substitute in Eqs. 11911 
and 121)11 for E the expression 12H and in the latter for a the expres- 
sion 1251 . I expand in powers of 5a and of the small entities £, {X' — X)'^ 
and (jijja) and keep terms up to the second order in Sa and up to first 
order in the other small entities. Some of the linear terms in 5a cancel 
because of the saddle-point condition. The exponent in Eq. 12011 takes 
the form 

aoL 





^ N 




^ Ne 


5a 


+ ^trg„ 


acSa 


+ 4a''"S. 



5 



'2A 



0-G(jl,j2 



+ ^(X-X')Vg, 



N 



acLSaL 
craScraGiiijh) 



Ne 



+-(x-x')Vg„ 



Ne 



aGSa{aGL) + aaSaaGLSaL 
Ne 



-trg,. 



+ ^(X-X')Vg, 



2A 

((7Gl/)^(o-G5cr)^ 



0"G(jl,j2)(o-G5cr)" 



(26) 



The leading terms in (Sa)^ are the first two terms in expression 12611 . 
These terms show that all massive modes have mass A''. We recall that 
N{X — X')^ is of order unity. Therefore, the remaining terms which are 
quadratic in 5a are negligible. The terms linear in Sa are all at most of 
order unity. To be non-negligible, they ought to be of order y/N. In the 
limit N ^ oo we are, thus, left with the first five terms in expression (I26II . 
This shows that the massive modes decouple from the Goldstone mode. 
Moreover, the massive-mode contribution attains exactly the form given 
in Ref. |^ and can, therefore, be integrated out without any problem. 
Hence, in contrast to the disorder problem studied in Ref. 0, the massive 
modes do not contribute to the strength of the parametric level correlation 
function in random-matrix theory. 
The result for Z is 



Z{E^,E2\X,X' ,.J) = 4 / d[cr]exp 



[lie , ^ , 
' 1 



N . 



+ -[X-X'YUg^[aGLy 



N^ 

\ trgc[(jij2)(TG_ 



8A2 



(27) 



I carry out the differentiation with respect to ji and j2 ■ The result is 
kie,X-X') = (l/2)y"d[a]exp| + ^trg„(aGL) 

+ ^(X-X')Vg„(aGL)^ 



x^(trgJ/(l)(aG)i.i])(trgJJ(2)(aG)2.2]) . (28) 



5 Integration over the Goldstone Mode 

Our result Eq. 1281 differs from the standard expression for the GOE two- 
point function by an additional term appearing in the exponent. Using 
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Eq. we rewrite this term in the form 

+ ^(X ^ X'ftrg^iaaLf = +^trgJ{[aG,L]f] . (29) 

Once again, the right-hand side of this equation shows very clearly that 
the term is due to the symmetry breaking caused by the perturbation. 

In the three graded traces appearing on the right-hand side of Eq. (I28II . 
the only matrices which break the pseudounitary symmetry are /(I) and 
/(2). Therefore, the two graded traces in the exponent depend only upon 
the "eigenvalues" (remaining integration variables). I use the parametriza- 
tions of both Refs. [5| and |3| to work out Z in the middle of the spectrum 
where Ao = 1. 

For the parametrization of Ref. I find 



k{e,X-X') oc / dAi / dXi J dX 

(l-A)A|Ai-A2| 



((1 + Ai)Ai(l + A2)A2)i/2(A + Ai)2(A + Aa)^ 
X exp<^ — ^(Ai + Az + 2A) 
7rr^ 

- — (Ai + Aa + 2A)(1 + Ai + Aa + 2A) 

x(Ai +A2 + 2A)2 . (30) 

For the parametrization of Ref. f4^, the two terms in the exponent take 
the form 

^(A - A1A2) - ^(2A?A^ - A? - A^ - A^ + 1) . (31) 

The first term agrees with Efetov's Eq. (5.35) if the definition of x follow- 
ing this equation is taken into account. The combination of integration 
variables appearing in the second term is the same as given by Simons 
and Altshuler. 

Hopefully, the derivation given above shows very clearly the role of 
symmetry breaking in parametric level correlations. The result confirms 
our expectation: The form of the parametric level correlation function is 
the same as for the disordered case. The strength factor differs and is 
given by 7rr^/(4d). Except for the numerical factor 7r/4, this result, too, 
corresponds to naive expectations: F^/d is a measure of the number of 
levels which are strongly mixed with each other as the external parameter 
changes from X to X' . In applications, this paramter can be estimated 
in terms of the strength of the perturbation and of the local mean level 
spacing. 



6 General Aspects of Symmetry Break- 
ing 

I now address more fully the symmetry-breaking mechanism which occurs 
when one considers parametric level correlations between two Hamiltonian 



7 



ensembles Hi and H2 (symbolically denoted by Hi < > H2). I do so in 

several situations with the intention of exhibiting the underlying similar- 
ities and differences. I consider the following cases: (i) GOE < > GOE; 

(ii) GUE < — > GUE; (iii) GOE < — > GUE. In the last case, the two Hamil- 
tonians Hi and H2 obviously do not belong to the same symmetry class. 
For the sake of comparison, I consider also (iv) the two-point autocorrela- 
tion function for the GOE — > GUE transition Proceeding as before, I 
calculate the resulting contributions to the effective Lagrangean (i.e., the 
additional symmetry-breaking terms in the exponent). These are jointly 
denoted by S and referred to as the parametric correlator. The commuting 
(c) and anticommuting (a) integration variables are arranged as follows: 
For the GUE two-point function, the sequence is (c, a, c, a) while for the 
GOE, it is (c, c, a, a, c, c, a, a). I define the graded matrix 

T3 = diag(+l, -1, +1, -1; +1,-1, +1,-1) . (32) 

Case (i) has been considered above. The parametric correlator was found 
to have the form S(i) = tTg{{[aG, L])^}. Case (ii) is formally very sim- 
ilar and leads to the same expression except that now aa and L have 
dimension four rather than eight. In case (iii), we get the GUE by adding 
to the GOE matrix in the advanced Green's function an imaginary ran- 
dom matrix. Thus, the new term in the exponent arises only from the 
advanced Green's function (and not from both, the retarded and the ad- 
vanced Green's function as in the previous cases (i) and (ii)). Moreover, 
the new term carries the matrix ra as the signal for the breaking of GOE 
symmetry and suppression of the Cooperon mode. As a result, the rele- 
vant term has the form ^(iii) — trg{([(TG, (Ig — Ls)T[i])^}. Case (iv) leads 
to a symmetry-breaking term of the form S'(iv) — trg{([(TG, Ts])^}. This 
obviously differs from ^(iii) . I observe that all these paramertic correlators 
have the same form, trg{ ([ctg, Tx])^}, with Tx given by 

Ti = 1/8 , 

Tii = Li , 
Tiii = (Is - L8)T3 , 

Tiv = T3 . (33) 

In summary, we have shown that the parametric correlation functions 
in random-matrix theory have a very simple form. Each one is obtained 
from the standard two-point function for level correlations by adding in 
the exponent of the generating function a term. That term is given by 
the graded trace of the commutator of the saddle-point solution ac with 
the particular matrix that describes the symmetry breaking in the actual 
case of interest. Except for a numerical factor which is of order unity, the 
factor in front of the commutator is given by /d. 
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